Anomalous Buckling of Charged Rods 
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Unscreened electrostatic interactions exert a profound effect on the onset of the buckhng instabihty 
of a charged rod. When this interaction is unscreened, the threshold value of the compressional force 
needed to induce buckling is independent of rod length for sufficiently long rods. In the case of rods 
of intermediate length, the critical buckling force crosses over from the classic inverse-square length 
dependence to asymptotic length-independent form with increasing rod length. It is suggested that 
this effect might lead to the possibility of the "electromechanical" stiffening of nanotubes, which 
would allow relatively long segments of them to be used as atomic force probes. 
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When slender objects are subjected to external com- 
pressional elastic forces, they are susceptible to bending 
deformations. The onset of such a deformation is known 
generically as the buckling instability P, Q ■ Mechani- 
cal failure also occurs in elastic and viscoelastic objects, 
such as falling ropes , "pouring" filaments of viscoelas- 
tic fluids Q as they hit horizontal solid surfaces, and mi- 
grating geophysical fluids reaching their terminations Q . 
An elastic rod of a given material and length can resist 
compressional forces up to a so-called critical buckling 
force Fc that increases with the (effective) bending stiff- 
ness of the rod, and decreases with an inverse-square law 
with its length L {Fc ^ \/L'^). Because of this length 
dependence, longer fllaments possess much lower critical 
buckhng forces 

The limitation on rod length imposed by the buckling 
instability is a major structural issue in nano-scale me- 
chanics (AFM-tips, nanotubes and nanorods 0, etc.) 
as well as the micro- and macro-scale. However, this 
limitation results from the local nature of elasticity, and 
may in principle be overcome if long-ranged interactions 
also exert a stiffening influence. Here we study the me- 
chanical response of an elastic charged rod to external 
compressional forces, and in particular the onset of Eu- 
ler buckling instability, taking into account the nonlocal 
nature of electrostatic self-interactions. For a cylindrical 
charged rod of radius r and surface number charge den- 
sity a , we find that long-ranged electrostatics leads, in 
the limit of a long rod, to a non-vanishing critical buck- 
ling force 

FciL -> oo) = A — eV^r^, (1) 

in which Eq is the permittivity of free space, e is the elec- 
tron charge, and A = 0.1137056 is a universal numerical 
prefactor. In the case of rods with a finite length L, we 
find that the above result smoothly crosses over to a lo- 
cal 1/L^ dependence as L decreases. Crossover to this 




FIG. 1: The elastic charged rod is hinged at two ends and is 
subject to a compressional force F. 

dependence occurs when L is small enough that the ac- 
cumulated electrostatic self-interaction has not yet over- 
whelmed local elasticity. We also determine the shape 
of the rod at the onset of buckling, and show that the 
buckling rod becomes considerably flatter in the interior 
as a result of electrostatic self-repulsion. 

The elastic charged rod is considered to be inexten- 
sible, and its energy consists of two contributions. The 
first part results from the elastic bending energy i?b ~ 
fJ^dsH{s)^. This energy is controlled by the intrin- 
sic bending modulus K of the elastic rod and contains 
no electrostatic contributions. The curvature H{s) is as- 
sumed to be a function of the arclength parameter s. For 
a homogeneous elastic rod of circular cross section that is 
made of a material with a Young's modulus E, we have 
K = Ir^E 0. 

The second contribution to the energy arises from elec- 
trostatic interactions, which can be written as E^x — 

frf dsds' ,s , for a rod whose conformation is 

2 JO |r(s) — r(s') I ' 

represented by a space curve r(s). The electrostatic cou- 
pling constant is defined as T = e^/(47reoa^), a being the 
average separation between neighboring charges along 
the line. Considering the cylindrical geometry of the rod, 
one can express the linear number charge density 1/a in 
terms of the surface number charge density a through 
a = {2TTra-)^^ , which yields T = (7r/£o)e^o'^r^. 

Finally, to study the onset of buckling upon apply- 
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ing a compressive force F, we add 
F ds cos 9 (s) to the total energy. 

One can estimate the critical buckling force for a neu- 
tral rod using a simple force balance argument. Imagine 
that the rod is bent into an arc of a circle of radius R 
with a corresponding arc angle 9, so that L = R6. Then 
we can calculate the bending energy -Eb(^) = 5^ = 
and the end-to-end distance x{9) = 2i?sin| = 
Lsin|/(|) as a function of the bending angle 9. The 
elastic force that resists bending at the onset of such arc 
formation can be found as Fb = —dE\^{9) / dx{9)\g=o — 
12K/L^, which slightly overestimates the exact critical 
buckling force Fco — tt'^K/L^ 0| due to the artificial 
assumption of constant curvature. A similar argument 
can be used to qualitatively account for Eq. |^ in the 
case of a charged rod with negligible intrinsic rigidity. 
Imagine that charges of unit magnitude are placed along 
the rod in a regular pattern at a distance a from each 
other (see Fig. The electrostatic repulsive force that 
the first charge experiences can then be calculated as the 
sum of the contributions from all the neighbors, namely, 
Fei(l) = T(l + l/22 + l/32 + ---) =7r2T/6. If the charged 
rod is to undergo compressional failure, the external force 
has to be greater than this Coulomb repulsion, thus yield- 
ing the scaling form for the critical buckling force as 
reported in Eq. To obtain the correct numerical 

prefactor, however, one should look collective failure cor- 
responding to the lowest threshold critical force. 

A convenient expansion of the total energy can be car- 
ried out in terms of a suitable deformation field. For the 
mechanical response that we would like to consider, it 
proves sufficient to focus only on planar deformations, 
which are conveniently characterized via the angle 9{s) 
that the rod's local unit tangent vector makes with its 
unperturbed orientation. Expanding the total energy up 
to second order then yields 8] 
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where the electrostatic kernel is given as 
C{s,s')= I dsi I ds2-, ^ — 

Jo Js^ {S2-S1Y 

X {(S2 - Si) [e(s - Si) - e(s - S2)] 5{s - s') 

- [e(s - Si) - e(s - S2)] [e(s' - si) - e(s' - s^)]] , 

(3) 

with Q{s) the Heaviside step function. Note that H{s) — 
d9{s)/ds in this representation. 

The total energy is controlled by the spectrum of the 
following total energy bilinear operator 

IC{F; s, s') = (-K-^ - f\ S{s - s') + TC{s, s'). (4) 
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FIG. 2: The critical buckling force for a charged rod as a 
function of its length (log-log plot) . (a) The solid line corre- 
sponds to K = lx 10"^^ Nm^ and T = 1 nN. (b) The dashed 
line corresponds to K = 100 x 10~^^ Nm^ and T = 1 nN. (c) 
The dash-dotted line corresponds to K — 1 x 10^^^ Nm^ and 
T = 100 nN. 



While the eigenfunctions of this operator 7/'„(s)'s are the 
same as those for /C(0; s,s'), the eigenvalues have the 
form of £n — F, if £n's are the eigenvalues of /C(0; s, s'). 
The onset of Euler instability is then found when the low- 
est eigenvalue Eq — F passes through zero. For example, 
if we switch off electrostatics by setting T = 0, the lowest 
eigenvalue will be Ktt'^/L'^ — F that goes non-positive at 
the critical force Fco — tt'^K/L'^. In Ref. [9|, we have 
been able to calculate the exact spectrum of the electro- 
static kernel £(s,s'). The eigenvalues are found to be 
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where 7 = 0.577216 is the Euler constant, and ip{x) = 
d\nT(x)/dx is the digamma function. Using this exact 
result, we can determine the critical buckling force for an 
infinitely long charge rod (or equivalently a charged rod 
with negligible intrinsic stiffness) as reported in Eq. 
above, in which the numerical prefactor is calculated as 
A = Ao = 7 -0(1/2) + 3/2 ~ 0.1137056. The existence 
of a non-zero critical force in this limit is a manifestation 
of the long-ranged nature of electrostatic interactions. 

Exactly at the onset of instability, the lowest eigenfunc- 
tion (the ground state) ipo{s) has a vanishing eigenvalue, 
and thus provides a nonzero solution to the homogeneous 
Euler-Largange equation corresponding to Eq. ||2Jl. With 
the appropriate boundary conditions, this solution pro- 
vides the shape of the elastic charged rod at the onset of 
instability. We use a cosine basis expansion for the eigen- 
function as tpo{s) = X]n=i ^« (^T^) corresponds 
to the boundary condition of a rod with two hinged ends. 
The shape of the deformed rod can be found as 

u(s) ~ [ ds'tpois') = — - sin f — - 
Jo ^ n-TT ^ ■ 



V L 



(6) 



Note that in the absence of electrostatics we have uo{s) = 
CO sin (^) 0. 

We study the spectrum of the total energy operator 
/C(0; s, s') numerically [l]| . and use it to find the critical 
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FIG. 3: The rescaled critical buckling force for a charged rod 
as a function of the charging parameter Q = TL^ /K. Note 
that three distinct series of data (hollow circles, triangles, and 
dashed line), corresponding to the different curves in Fig. |5| 
have been collapsed on top of a universal curve. The solid 
line represents the interpolation formula of Eq. (|7J for com- 
parison. 



buckling force and the shape of a charged rod of arbitrary 
length at the onset of compressive failure. In Fig. |21 
the critical buckling force is plotted as a function of the 
length of the rod, for various values of K and T. The plot 
shows an small-L dependence for the critical force 

that crosses over to an L-independent saturation value 
as L passes through a crossover length scale i^, where 
cx/??7T. 



The critical buckling forces corresponding to different 
values of the parameters K, T, and L can be collapsed 
onto a universal curve as shown in Fig. |31 if normalized 
with the critical buckling force of the neutral chain Fco 
and plotted as a function of the dimensionless charging 
parameter Q = TL^ / K. An interpolation formula of the 
form 
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is found to satisfactorily represent the universal curve as 
revealed by the comparison in Fig. O 

The shape of the charged rod at the onset of buckling 
is also calculated, and shown in Fig. 2]for various values 
of the charging parameter Q. It appears that charging 
leads to deviations in the shape of the buckling rod from 
the sinus-profile |2jJ in that there is enhanced flattening in 
the interior. This is to be expected because the interior 
of the charged rod is subject to stronger build-up of elec- 
trostatic self-repulsion as compared to the end-segments 
where "half" of the repelling charges are absent. A simi- 
lar effect has also been observed in the bending response 
of charged elastic rods . 

We can also study the buckling instability for charged 
rodlike polymers, or polyelectrolytes, in a solution by us- 
ing the screened Debye-Hiickel interaction e~'^^ /{er) in- 
stead of the long-ranged Coulomb l/r interaction, where 
e is the zero-frequency dielectric constant of the solu- 
tion and represents the Debye screening length. The 




FIG. 4: The shape of a charged rod at the onset of Euler buck- 
ling instability. The dashed line corresponds to Q = 0, the 
solid line corresponds to Q = 10^, and the dash-dotted line 
corresponds to Q = 10®. The buckling charged rod flattens 
in the interior as the charging is increased. 



eigenvalues of the total energy can be calculated numeri- 
cally for the screened case, and can be used to determine 
Fc for each set of parameters. We find that for sufficiently 
long polyelectrolytes the critical buckling force is given as 
Fc = TT^kBT{£Q + £osf)/L^, in which £o = K/ksT is the 
intrinsic persistence length, and ^osf = T/{AekBTK^) is 
the well-known Odijk-Skolnick-Fixman electrostatic per- 
sistence length for polyelectrolytes d, 0|. This result 
confirms that the so-called wormlike chain description of 
polyelectrolytes remains valid in determining their me- 
chanical response and onset of buckling instability, as 
long as kL ^10, ^'^^ shorter chains, however, 

the dependence is altered, and a crossover similar 

to the one described in Fig. |21 takes place. As a par- 
ticular case of interest for assessing the wormlike chain 
model of polyelectrolytes one can consider the unscreened 
limit (k = 0), where it predicts a crossover of the form 
= 4^ + f^T = Fco (1 + ^if^Q) 0, which is to 



be contrasted with Eq. (UJ above. 

The familiar image of a long-haired girl touching the 
van de Graaff machine suggests that a practical way of 
imposing the required charging is by applying a voltage. 
For a conducting cylinder of length L and radius r that is 
kept at a potential V relative to "infinity" one can 
calculate the induced surface charge density, and deduce 
from it the corresponding asymptotic critical buckling 
force as 



F,(L^oo,y) = 



[HL/r)f 



(8) 



For a thread of human hair we have r ~ 0.1 mm and 
A\air - 10"" Nm^ which yields for L = 1 cm a critical 
force of Fc ~ 10"'' N. Applying a voltage of ^ = 50 kV 
(typical of van de Graaff generators) then results in a 
critical force of Fc ~ 10~^ N for a one-meter long piece 
of hair! 

Perhaps the most interesting venue in which these re- 
sults find application will be in hardening of atomic force 
probes. Carbon nanotubes have been found to be struc- 
turally quite robust and have exceptionally high Young's 
modulus (in the TPa range) [la|. However, the fact 
that they can grow up to microns in length while having 
nanometric diameters renders them quite susceptible to 
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buckling. The buckling of multi-walled carbon nanotubes 
have been recently investigated experimentally by Dong 
et al. [l^ . In their experiment, a 6.9 /zm long nan- 
otube have been placed under compression and its criti- 
cal buckling force (typically in the nN range for micron- 
sized nanotubes) have been measured, from which they 
could extract the bending rigidity of the nanotube as 
ifnanotubc = 8.641 X lO^^o Nm^ JJ,] . Using Eq. ©, we 
can now estimate that a carbon nanotube 30 nm in diam- 
eter can resist forces up to ^ 1 nN even when it is 1 mm 
long, provided it is kept at a voltage of 200 V. This im- 
plies a remarkable "electromechanical stiffness," in con- 
trast with the intrinsic mechanical resistance to buckling 
which is diminished by a factor of 10^ when the length 
of the nanotube is increased from a micron to a millime- 
ter. The stiffening mechanism may also be useful in the 
recently reported nanometer-scale electromechanical ac- 
tuator that is based on a multi-walled carbon nanotube 
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